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Abstract 

The investigation of the oscillation pattern induced by the sterile neutrinos 
might determine the oscillation parameters, and at the same time, allow to 
probe CPT symmetry in the leptonic sector through neutrino-antineutrino 
mass inequality. We propose to use a large scintillation detector like JUNO 
or LENA to detect electron neutrinos and electron antineutrinos from MCi 
electron capture or beta decay sources. Our calculations indicate that such 
an experiment is realistic and could be performed in parallel to the current 
research plans for JUNO and RENO. Requiring at least 5cr conhdence level 
and assuming the values of the oscillation parameters indicated by the current 
global £t, we would be able to detect neutrino-antineutrino mass inequality 
of the order of 0.5% or larger, which would imply a signal of CPT anomalies. 

Keywords: Sterile neutrino. Neutrino oscillometry, CPT, Liquid 
scintillator detector. Artificial source 


1. Introduction 

The discovery of neutrino oscillations [1] provided undisputed evidence of 
incompleteness of the conventional Standard Model (SM) of particle physics 
[2] and gave a clear indication that further neutrino studies may lead us 
towards phenomena Beyond the Standard Model (BSM). One of such possible 
manifestations would be violation of CPT invariance [3] in the process of 
neutrino oscillations. 
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CPT symmetry is considered a fundamental law of nature, to which we 
know no exceptions yet. One of the important implications of this symmetry 
law is the requirement that in all processes involving particles and their 
antiparticles the corresponding probabilities and other parameters governing 
the process should be exactly the same. Many experiments based on that 
principle have been conducted to probe the CPT conservation. For instance, 
in the leptonic sector, a measurement of ^'-factor for electron and positron 
performed at the University of Washington in Seattle jl] or, in the baryonic 
sector, a mass measurement of proton and antiproton implemented at CERN 

[SI El- 

Neutrinos carry important information about weak interactions. Espe¬ 
cially their interferometric nature, make them a very sensitive tool to search 
for unconventional physics, in particular CPT violation and Lorentz violation 
jam. However, until now there was no direct and simultaneous experiment 
with comparison of neutrino and antineutrino, therefore the issue of CPT vi¬ 
olation in the leptonic sector remains unresolved. The mass inequality could 
give a unique possibility to do it. 

Obviously, the properties of neutrinos doom the attempts to verify CPT 
conservation through oscillation measurements to a likely failure or, at best, 
have a prohibitively high price tag. Very low interaction cross-sections, lim¬ 
ited neutrino flux, and the size of the available detectors exclude such at¬ 
tempts on short time scales. However, if the active neutrinos like the well- 
known electron-, muon-, and tau neutrino and their antiparticles were not 
the only existing neutrino flavours, it might be possible to search for CPT 
violation already in the near future, as it is outlined in this article. 

One of the exciting new hypotheses is the existence of the light sterile 
neutrinos with unknown nature and properties jS]. The existence of one or 
more sterile flavors is indicated by the re-evaluation of reactor neutrino fluxes 
experiments [9], some of the beam experiments such as MiniBoone, LSND 
and the Gallium anomaly mm- 

Although the direct observation of sterile neutrinos is not possible, their 
presence would be reviled by a distinct oscillation pattern, caused by the con¬ 
version from active to sterile flavor. If observed, the pattern would also give a 
hint of the number of sterile neutrino flavors. The simpliest assumption is the 
scenario with one additional sterile flavour, so called 3-1-1 scheme m- The 
current global £t to the oscillation data suggests that the additional mass- 
squared difference and mixing angle governing the active-to-sterile conversion 
should be of the order of Amh 1 eV^ and sin^ 29ee 0.1 respectively n- 
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Assuming these parameter values and relying on the proportionality between 
the energy and oscillation length, we may choose the energy of neutrinos in 
a way that at least one full oscillation period will overlap with the active vol¬ 
ume of a single detector, like it is proposed in the Borexino SOX experiment 
|13j . The realization of this neutrino oscillometry method [H] will open a 
unique opportunity to search not only for disapperance oscillations to sterile 
neutrinos but also to check the validity of CPT conservation in the leptonic 
sector if such a transition is observed. 

2. Probing neutrino-antineutrino mass inequality 

If CPT is strictly conserved the oscillation of neutrinos and antineutrinos 
would be described by identical probability functions with two main param¬ 
eters: mixing angles and mass-squared differences. For instance, in the case 
of 3-1-1 scenario, the survival probability of electron neutrino and electron 
antineutrino can be approximated with [T^ : 

p(v\ -)■ ^v\) = 1 - sin^(26'ee) ■ sin^(1.27 ■ ■ ^). (1) 

tj 

Note that the equation ([^ is valid for short baseline experiment (SBL). 
Consequently, if CPT conservation holds, we should obtain in the classical 
case [T5] : 

(Amli)jv ^ 

(Amli)^ ’ ^ ^ 

where N stands for the value obtained from the measurement with neutrinos 
and A, with antineutrinos. At the same time, if in the real experiment 
the mass inequality will be observed, then it would mean CPT violation 
and the estimation of this process should be done within the framework of 
the Standard Model Extension (SME) formalism with using coefficients of 
Lorentz- and CPT-violation [3]. In this paper, the comparison is made only 
for the ratio of the mass-squared differences, because Am^^ can be extracted 
from the experimental data with better precision than ©ee. For instance, 
it depends far less on the intensity normalization. The evaluation of ©ee 
will be considerably less accurate, because of the systematic uncertainties in 
absolute values of neutrino flux and hducial volume. 

It should be noted, that P{i'e —t z/e) ^ P{i^e —t would mean both 
CPT and CP violation. However the latter can not be probed in survival 
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experiments if CPT is conserved]^ 

3. Outline of the experiment 

Successful application of neutrino oscillometry requires some degree of 
energy selection and well-defined source position. For these reasons the ex¬ 
periment cannot rely on natural neutrino sources such as the Sun, the Earth 
or cosmic neutrinos. Nuclear power reactors also are excluded, as they do not 
produce considerable amounts of electron neutrinos. The best alternatives 
are provided by man-made, high-intensity beta decaying radioactive sources 
emitting electron neutrinos and antineutrinos with the energy around 1 MeV. 
This choice is based on the outcome of the global £t predicting of the 

order of 1 eV^. Precise measurement of that value both for electron neutrinos 
and antineutrinos will be the key outcome of the proposed experiment. 

The mass-squared difference will be derived from the reconstructed os¬ 
cillation probability curves. To measure the oscillation pattern the detector 
has to have sufficiently low energy threshold for the registration of both neu¬ 
trinos and antineutrinos and have adequate position and energy resolution. 
Furthermore, to achieve the required statistics, the active volume of the de¬ 
tector should ideally extend over at least one oscillation period. The only 
devices that can currently fulhl these requirements are large liquid scintillator 
detectors. 

3.1. Neutrino and antineutrino sources 

The obvious hrst choice for a z/g source is ^^Cr, well known from GALLEX 
Cain] and GNO [TH] experiments. The decay scheme of ®^Gr is shown in 
Fig.ll This nuclide has 27.7 days half-life and a neutrino spectrum dom¬ 
inated by mono-energetic lines at an energy of around ~0.75 MeV (90.1% 
total branching ratio), as shown in Fig. Gurrent estimations for the SOX 
project [13] indicate that a 35 kg chromium oxide source enriched in ^‘^Gr 
(~30%) placed in a nuclear reactor, can be activated up to 10 MGi. Ac¬ 
counting for the decay during extraction from the reactor, transport, and 
installation inside the detector, the realistic activity of the source at the 
start of the experiment would be around 8 MGi. 

As the source of Pf. we propose ^^^Ge —Pr chain [20]. The decay scheme 
is shown in Fig. The antineutrino energy spectrum of ^^^Pr is continuous 


^suggested by E.Kh. Akhmedov 
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(27.7 days) 



Figure 1. The decay scheme of ®^Cr, 
with indicated energies of emitted neu¬ 
trinos |22) . 



Figure 2. The decay scheme of system 

i44cg_i44 Pr [1^12^. 


with the end point energy at about 3 MeV [21] and with the overall half-life 
of 285 days. About 48.5% of the emitted antineutrinos are at energies above 
the detection threshold for the IBD reaction (see Section 3.2) and can be 
used in the experiment. Current studies indicate that an activity of 100 kCi 
or more can be reached for the Ce-Pr source [23] . In our calculations we 
assume the initial activity of 120 kCi. 

The properties of both sources are summarized in Table 


Table 1: Properties of sources proposed as neutrino- and antineutrino emitters for the 
oscillometry measurement. 


Source 

Type 

Spectrum 

E, MeV 

A, MCi 

Ti/ 2, d 

Expos., d 

51 Qr 


mono 

0.75 

8 

27.7 

55 

144Qe -144 pr 

Ve 

continuous 

1.8 - 3.0 

0.12 

285 

300 


3.2. Neutrino detection in a liquid scintillation detector 

Elastic neutrino-electron scattering is the dominant interaction channel 
leading to the registration of Ug events in a liquid scintillator detector. This 
reaction 

Z/e + e“ —)■ 
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has similar characteristic to the Compton effect. While all three neutrino 
flavours interact with electrons by neutral current, electron neutrinos interact 
as well via charged current. The cross-section of this interaction may be 
approximated by a simple expression: a = 0.7 x Ej, ■ 10“^'^cm^ [2l]. The 
uncertainty of this estimation, including of radiation corrections, is about 
2 % 121 ]. 

The dominant detection channel of electron antineutrinos in liquid scin¬ 
tillator is the inverse beta decay reaction (IBD) 

Pe + P —^ ^ 

This reaction produces two distinct signals in the detector allowing for an 
unambiguous identihcation of a neutrino event. The kinetic energy and anni¬ 
hilation of the positron and the second, from neutron capture on hydrogen. 
The former is prompt and contains the information on the energy of the 
neutrino. The later 2.2 MeV gamma signal is delayed by about 250 /is. The 
threshold energy for IBD is 1.8 MeV and the cross section of the process can 
be approximated by the expression a = 9.5 x {E^ — 1.29)^ ■ 10“^^cm^ [25] . 

3.3. Proposed implementation and expected event rates 

To evaluate the feasibility of the proposed experiment the event rates were 
calculated for two alternative shapes of liquid scintillator detector: spheri¬ 
cal, with the source installed in the centre, and cylindrical, with the source 
mounted close to the center of the end cap of the tank. 

The specihcations of the spherical detector are those of JUNO [2S], cur¬ 
rently under construction in China. It will contain 20 ktons (hducial) of 
liquid scintillator enclosed in a spherical volume with a diameter of 34.5 m. 
Spherical shape provides the optimal detection geometry when the source of 
z/e or Pe is placed in the center, as shown on Fig. j^. The expected exposition 
time is restricted by the half-life. For the Cr source we therefore assume 55 
days and for the Ce-Pr source 300 days. Gamma-ray background from the 
source inside the detector can be reduced by using a radiation shielding, e.g. 
from tungsten [T3|. 

To evaluate the case of a cylindrical detector, we have used LENA speci¬ 
fications [22] : 100 m high, 28 m diameter, 50 kton hducial mass. The source 
of z/e or Pe would be located close to the top of the detector, as it is shown on 
Fig. [^. In this case only the hemispherical volume near the top would be 
considered for analysis. This is clearly a less favourable geometry as it uses 
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A. B. 

Figure 3. The scheme of the proposed experiment. Left panel: the configuration for JUNO 
with the source in the center of the detector. Right panel: the configuration for LENA 
with the source in the center of the top of the tank. 


only part of the fiducial volume and has to cope with non-optimal coverage 
by the photo sensors (PMTs). The only considerable advantage of this ge¬ 
ometry is that it eliminates the need to immerse the source into the volume 
of the detector. For the purpose of this analysis the two detector geometries 
differ just by a factor of 4 (volume of a sphere vs. hemisphere) hence energy- 
and distance-dependent event rate in both cases can be expressed by the 
equation |20l EH] : 

Emax 

/ 

Emin 

where Aq is the source activity at the start of the measurement, n the density 
of free protons (electrons) in the target, A the width of bin with the center 
L, (t{E) the cross section for IBD or electron-neutrino scattering, S{E) the 
spectral shape (delta-function for monoenergetic neutrinos), P{L,E) the os¬ 
cillation probability ([^, te the time of measurement and A the decay constant 
of the source. 

The total exposure time needed to complete the measurement would be 
about 355 days assuming that both sources are used subsequently. The event 


N{L,E)= j 


^■n-a{E)-S{E)-P{L,E)-{l- 

A 


■exp[—Ate]) di? dL, (3) 
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rates were plotted as a function of the distance from the source and are shown 
on Fig. 1^ and Fig. for i/g and z/g respectively. The presented calculations 
correspond to the spherical geometry. 




A. B. 

Figure 4. The expected rate in JUNO detector. Left panel: the total rate for ^^Cr source, 
time of measurement 55 days. Right panel: the total rate for Pr source, time 

of measurement 300 days. 

Since the detection channels for z/g and z/g are different and disentangled, it 
would be possible to make the measurement with both radioactive sources at 
the same time.This conhguration would even have the advantage of reduced 
systematic error due to normalization. However, from the implementation 
point of view, that would be a major experimental challenge. 


4. Results 


The main outcome of the experiment will be the independent determina¬ 
tion of the oscillation parameters for neutrinos and antineutrinos by htting 
the spectra shown on Fig. As it was explained earlier, the main parameter 
for testing CPT symmetry is Amfi- The conhdence level can be evaluated 
from equation ([^, 


(Ami^)jv ± aN 

{Aml^)A ± (Ta 


= 1 ± CT, 


(4) 


where ctat is the uncertainty of measurement of {Amh)N, cta the uncertainty 
of measurement of a the total uncertainty of ratio. Since the 
























starting assumption is that {/S.m\^)]si and are identical (although 

with different error bars), the relative uncertainty can be estimated as follows: 

- {Ami){Ami) a' ^ ’ 

The event samples were generated by using Monte-Carlo approach. The 
basic steps which were used in analysis: 

• The oscillation spectrum was generated with certain oscillation param¬ 
eters {Ami ^^cl sin^26*ee)- 

• The pattern of real oscillation spectrum was obtained with adding of 
basic uncertainties: statistical error, the error of position resolution, 
the error of source activity and the error of fiducial volume of detector. 
For antineutrino case the energy resolution was also used in analysis. 

• The oscillation parameters were extracted from the pattern of real spec¬ 
trum with using fit-function the same as probability function 0 taking 
energy and position resolution smearing into account. 

• The £t-function was optimized by using the method of minimization 

from package Minuit2 ROOT. 

• The gaussian distribution for uncertainties of oscillation parameters 
was assumed. 

The example of such distribution of errors is shown on Fig. As can be seen, 
the relative error for sin^ 29ee is larger at least by one order of magnitude than 
the relative error for Ami- It explains also our choice for Ami 
paprameter parameter of OPT test. The uncertainty of these measured values 
of Ami determined independently for both sources, and combined with 
the equation ([^. However this approach allows also to determine ©ee with 
high precision. 

For detector-related uncertainties, the design parameters given by JUNO 
and LENA collaborations were used. In the case of JUNO the energy reso¬ 
lution of 3%/ E[MeV] and position resolution of 9 cm.{vl) and 4.5 cm(z/e) 
were used. For LENA 6.1%/ E[MeV], 10 cm(i/e) and 4.5 cm(Pe) were used 
respectively. 

Neglecting background from radioactive decays, the main background in 
the energy region of interest for Ve detection arises from the solar neutrinos 
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sin''20„ = O.1 


Ami, - ^ 




measured value - true value 


A. 


B. 


Figure 5. The distribution of difference between reconstructed value and true value for 
neutrino case in LENA detector geometry. Left panel: the reconstructed parameter is 
Amh with true value 1 eV^. Right panel: the reconstructed parameter is sin^ 20ee with 
true value 0.1. 


(^Be). It is shown in Table The inflnence of this factor on the hnal resnlt 
was also taken into acconnt in onr estimations. Nevertheless, direct mea- 
snrements of the solar backgronnd are carried ont antomatically whenever 
the sonrces are removed so the correction may be done afterwards. In the 
case of antinentrinos, the highest backgronnd comes from nnclear reactors 
and can be neglected, even in the case of JUNO detector, becanse the ex¬ 
pected nentrino flnx from ^^^Ce —Pr sonrce located in the middle of the 
detector is considerably higher than that of 12 power reactors 53 km away. 
The estimated event rates from sonrces with the strength of 8 MCi for ^^Cr 
and 0.12 MCi for ^"^^Ce Pr as well as the corresponding backgronnds are 
presented in Table 

Table 2: Expected background and signal events for full time of the experiment. The mass 
of JUNO’s sphere is 20 kt, mass of semisphere in LENA is 5 kt. 


Type 

Signal 

JUNO 

Backgronnd 

JUNO 

Signal 

LENA 

Background 

LENA 

Exposition 
time, d 


1148860 

348032 

459234 

87008 

55 


466800 

3090 

193500 

45 

300 


The main ontcome of onr simnlations is shown in Fig. (for LENA) and 
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Fig. 0 (for JUNO). The figures show a 5a sensitivity limit calculated for 3 
different sin^ 26^66 values and plotted as a function of Amji. For instance, if 
we check the performance for the parameters chosen by the global £t that 
is > 1 eV^ and sin^26'ee ~ 0.1, we conclude from Fig. that LENA 

will detect (at the 5a level of measured events) a difference of neu- 

trino/antineutrino and therefore CPT violation that exceeds 1%. For the 
same parameters JUNO would reach 0.5% sensitivity as it may be observed 
from Fig. The slightly better sensitivity of the JUNO detector is a conse¬ 
quence of the spherical shape and better expected energy resolution. 



Am^i [eVj 


Figure 6. Sensitivity limits of LENA for the detection of mass inequality at 5a level of 
measured events calculated for 3 different mixing angles and plotted as a function of Amh. 
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Figure 7. Sensitivity limits of JUNO for the detection of mass inequality at 5cr level of 
measured events calculated for 3 different mixing angles and plotted as a function of . 
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5. Conclusions 

We have outlined a scheme to probe the mass inequality by analyzing the 
neutrino and antineutrino oscillation patterns induced by the still hypothet¬ 
ical sterile neutrinos. If the ’’sterile” oscillation patterns are detected, the 
experiment will yield new fundamentally important oscillation parameters 
and will be sensitive to the mass inequality on the level of a few tenths of 
a percent at 5a signihcance. This sensitivity can be achieved assuming the 
current best-£t values of the oscillation parameters indicated by the global 
£t. 

The proposed experiment would utilize a large scintillator like future 
JUNO, LENA or RENO-50 experiments. Since JUNO is now approved and 
has just entered the hrst construction phase, the proposal is certainly feasible. 
The same concerns the use of MCi sources. Since similar ones were already 
produced and used by GALLEX and GNO experiments, there is already the 
necessary knowhow and, in the case of Gr source, also enriched material in 
sufhcient quantity. 
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